FRICTION  ON  CRACK  SURFACES  DURING  COMPRESSION  OF  EXPLOSIVES  -  A 
POSSIBLE  IGNITION  SOURCE  FOR  UNPLANNED  EXPLOSIONS  DUE  TO  ACCIDENTAL 

DEFORMATIONS 


Donald  A.  Wiegand*,  Brett  Redingius 

•  ARDEC  Picatinny,  NJ,  USA. 

•  Kevin  Ellis  and  Claire  Leppard 

Atomic  Weapons  Establishment,  Aldermaston,  UK 


ABSTRACT 

The  mechanical  properties  of  composite  plastic 
bonded  explosives  have  been  investigated  as  a  function 
of  confining  pressure.  The  results  indicate  different 
failure  processes  in  two  pressure  ranges,  a  low  pressure 
range  between  about  0.1  to  7.0  MPa  which  is  considered 
in  this  paper  and  a  higher  pressure  range.  In  the  low 
pressure  range  crack  processes  are  important  in  failure. 
The  pressure  dependence  of  the  compressive  strength  in 
the  low  pressure  range  is  attributed  to  coulomb  friction 
between  surfaces  of  closed  shear  cracks  and  from  the 
observed  linear  increase  of  the  strength  with  pressure  and 
the  angle  of  the  fracture  plane  a  friction  coefficient  is 
obtained.  A  friction  coefficient  can  also  be  obtained 
from  the  ratio  of  the  compressive  to  tensile  strength  and 
directly  form  the  above  angle.  The  friction  coefficients 
obtained  from  these  three  separate  observations  are  in 
agreement  and  this  is  taken  as  strong  evidence  for  the 
importance  of  this  friction  in  determining  strength  and 
mechanical;  failure.  Frictional  heating  during 
deformation  can  then  cause  hot  spots  leading  to  ignition  . 


1.  INTRODUCTION 

Explosives  and  propellants  are  often  used  under 
conditions  of  confinement  and  pressurization. 
Explosives  are  confined  in  projectile  cases  and  are 
pressurized  during  launch  by  set  back  forces  and  during 
impact  by  set  forward  forces.  Propellants  are  confined 
by  the  breech  and  are  pressurized  by  hot  gases  during 
burning.  Because  of  these  pressurizations  the  properties 
of  explosives  and  propellants  under  pressure  are  of 
interest.  In  particular  the  mechanical  properties  under 
pressure  are  needed  for  modeling  and  safety 
considerations.  For  example,  for  the  modeling  of  an 
explosive  filled  projectile  during  launch  or  impact  the 
mechanical  properties  of  the  explosive  under  pressure  are 
required.  Similarly,  for  the  modeling  of  a  propellant 
charge  during  burning  the  mechanical  properties  of  the 
propellant  under  pressure  are  required.  Of  particular 
concern  are  the  mechanical  failure  properties  under 
pressure.  Fracture  or  yield  during  the  use  of  explosives 
can  lead  to  unwanted  and/or  hazardous  ignitions  (Howe 
et  al,  1985)  (Frey,  1985)  (Coffee,  1985)  (Heavens  and 


Fields,  1972).  In  addition,  the  fracture  of  propellants 
during  burning  can  lead  to  hazardous  burning  conditions 
(Nicolaides  et  al,  2000).  The  results  presented  here  also 
indicate  the  possible  hazards  associated  with  crack 
processes  (Wiegand  et  al,  2008).  Because  of  these 
considerations  a  program  has  been  initiated  to  study  the 
mechanical  properties  of  these  materials  under  hydrostatic 
pressure  (Wiegand,  2000a,  2000b)  (Wiegand  and 

Reddingius,  2003,  2005a,  2005b). 


2.  EXPERIMENTAL 

A  high  pressure  chamber  designed  at  Structural 
Behavior  Engineering  Laboratory,  Phoenix,  AR  to  contain 
pressures  up  to  138  MPa  was  used  to  study  the 
compressive  mechanical  properties  as  a  function  of 
confining  pressure  (Wiegand,  2000a).  Hydraulic  oil  was 
used  as  the  confining  medium  and  the  sample  in  the  form 
of  a  right  circular  cylinder  was  protected  from  the  oil  by  a 
tight  fitting  tubular  gum  rubber  or  neoprene  shroud.  The 
ends  of  the  sample  were  against  steel  platens  and  O-ring 
seals  were  used  to  prevent  oil  from  reaching  the 
sample. (Wiegand,  2000a,  2000b)  (Wiegand  and 

Reddingius, 2003,  2005)  The  confining  pressure  is  taken 
here  as  the  chamber  hydrostatic  pressure  before  the  start  of 
and/or  during  the  axial  compression.  In  all  cases  the 
pressures  referred  to  here  are  this  hydrostatic  pressure. 
The  chamber  pressure  was  determined  using  a 
SENSOTEC  pressure  gauge,  model  JTE/1 108-03, 
calibrated  by  the  manufacturer  and  mounted  at  the  base  of 
the  chamber.  In  addition,  A  McDaniel  Controls  dial 
pressure  gauge  was  mounted  at  the  pump.  Measurements 
at  atmospheric  pressure  (0.1  MPa)  were  made  in  air. 

The  samples  were  compressed  along  the  cylindrical 
axis  and  two  LVDT's  (linear  voltage  differential 
transformers)  were  mounted  to  measure  axial  strains. 
They  were  spaced  180  degrees  apart  around  the 
circumference  of  the  sample  with  their  axes  parallel  to  the 
sample  axis.  The  sample  axial  strain  was  taken  as  the 
average  of  the  strains  obtained  from  the  two  LVDT's.. 
Two  or  three  additional  LVDT's  were  mounted  to  measure 
radial  strains.  They  were  placed  in  a  plane  at  the  sample 
axial  mid  position  with  their  axes  perpendicular  to  the 
sample  axis.  They  were  also  180  degrees  apart  (or  120 
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TABLE  1 


Composition  of  Composites  Discussed 


Binder 


Name 

Explosive/Inert 

Polymer 

Plastizer 

TG(°C) 

EDC37 

HMX 

NC 

DNEB/TNEB 

-63 

91% 

1.0% 

5.22%/2.78% 

PAX2A 

HMX 

CAB 

BDNPA/F 

-37 

85% 

6% 

9% 

PBS  9501 

SUCROSE 

ESTANE 

BDNPA/F 

-41(B) a 

94% 

3% 

3% 

PBX  9501 

HMX 

ESTANE 

BDNPA/F 

-41(B)a 

95% 

2.5% 

2.5% 

LX-14 

HMX 

ESTANE 

-31(B)a 

95.5% 

4.5% 

Nomenclature:  HMX  -  Cyclotetramethylene  tetranitramine..  NC  -  Nitrocellulose.  DNEB  -  Dinitroethylbenzene,  TNEB  - 
Trinitroethylbenzene.  CAB  -  Cellulose  Acetate  Butyrate.  BDNPA/F  -Bis(2,2-Dinitropropyl)Acetal/Formal.  Estane  -  Polyurethane. 
B  -  Property  of  the  Binder, 
a  (Wiegand  and  Reddingius,  2005) 


degrees  for  three  radial  LVDT's)  around  the  sample 
circumference  (Wiegand,  2000a,  2000b)  (Wiegand  and 
Reddingius, 2003,  2005). 


Axial  stress  versus  axial  strain  data  in  compression 
were  obtained  using  the  above  chamber  and  an  MTS 
servo-hydraulic  system  operated  at  a  constant 
displacement  rate  (Wiegand,  et  al,  1991)  (Pinto  et  al, 
1985).  All  work  presented  was  carried  out  at  strain  rates 
of  approximately  0.0005/sec.  and  0.001/sec  with  the 
exception  of  results  explicitly  a  function  of  strain  rate. 
The  right  circular  cylinder  samples  were  3.81  cm  (1.50 
inch)  in  length  and  1.90  cm  (0.75  inch)  in  diameter  and 
so  had  a  length  to  diameter  ratio  of  two.  The  end  faces 
of  all  samples  were  coated  with  a  lubricant  to  minimize 
frictional  effects  between  the  sample  end  faces  and  the 
loading  platens.  The  sample  temperatures  during 
measurements  were  between  20  and  23  C  and  samples 
were  conditioned  at  temperature  for  at  least  two  hours 
before  measurement.  The  dimensions  of  all  samples  at 


0.1  MPa  (atmospheric  pressure)  were  used  to  obtain 
engineering  stress  and  engineering  strain. 

Most  of  the  measurements  reported  here  were  made 
with  samples  of  EDC37,  a  UK  plastic  bonded  explosive 
(PBX).  The  composition  of  this  explosive  is  given  is 
Table  1  along  with  the  compositions  of  other  similar 
materials  which  are  referred  to  in  this  paper  (Wiegand  and 
Reddingius).  Samples  were  prepared  by  pressing  into 
large  billets  and  machining  to  size.  Precautions  were  taken 
to  insure  that  the  cylinder  end  faces  were  adequately  flat 
and  parallel  (Wiegand,  1998).].  The  densities  of  all 
samples  were  in  a  narrow  range  close  to  the  maximum 
theoretical  (zero  porosity)  density.  The  densities  of  most 
samples  were  determined  before  and  after  compression  by 
weighing  in  air  and  in  purified  water  and  by  using  the 
density  of  water  at  the  temperature  of  measurement.  All 
sensors  were  calibrated  by  the  manufacturer  or  calibrated 
against  standards  provided  by  the  manufacturer.  It  is 
estimated  that  variations  from  sample  to  sample  in  any 
measured  quantity  are  significantly  greater  than  errors 


introduced  by  the  sensors  or  errors  introduced  during 
data  processing. 


3.  RESULTS 

The  results  indicate  two  pressure  ranges  in  which  the 
mechanical  failure  properties  differ,  a  low  pressure  range 
between  about  0.1  and  7.0  MPa  for  EDC37  in  which 
failure  occurs  via  crack  processes  and  a  higher  pressure 
range  between  about  7.0  and  138  MPa  for  the  same 
composite,  in  which  failure  occurs  via  slip  processes. 
The  low  pressure  range  is  here  considered  in  some  detail. 
The  higher  pressure  range  is  considered  elsewhere. 
(Wiegand,  et  al,  to  be  published). 

In  Figure  1  axial  stress  versus  axial  strain  curves  for 
EDC37  are  given  for  several  confining  pressures  in  the 
low  pressure  range.  The  pressures  are  as  marked.  The 
curve  at  atmospheric  pressure  (0.1  MPa)  has  the  typical 
response  of  many  energetic  materials  at  this  pressure, 
and  temperature,  i.e.,  an  initial  linear  increase  of  stress,  a 
maximum  stress  followed  by  work  softening  (decrease  of 
stress)  with  increasing  strain  (Wiegand,  1985,  1998). 
Available  evidence  indicates  that  the  deviation  from 
linearity  after  the  initial  linear  region,  the  maximum  and 
the  work  softening  are  associated  with  crack  damage 
(Wiegand,  1985,  1998).  In  fact  acoustic  emission  results 
suggest  that  some  crack  damage  is  generated  throughout 
the  stress-strain  curve  (Ellis  et  al,  2005).  Evidence  of 
surface  cracking  is  observed  for  most  samples  deformed 
into  the  higher  strain  parts  of  the  work  softening  region 
at  atmospheric  pressure. 

Both  the  initial  slope  and  the  maximum  stress 
increase  with  increasing  pressure  while  the  negative 
work  softening  slope  decreases  with  increasing  pressure 
and  is  close  to  zero  at  a  pressure  of  6.9  MPa.  The  latter 


EDC37 


Figure  1.  Axial  stress  versus  axial  strain  for  pressures  in 
the  low  pressure  range. 


result  suggests  that  either  crack  damage  decreases  with 
increasing  pressure  and  is  negligible  at  6.9  MPa  or  that 
crack  damage  is  generated  but  that  other  factors  prevent 
the  work  softening  from  being  observed.  The  latter  is 
apparently  the  case  as  discussed  elsewhere.  (Wiegand  et  al, 
unpublished)  Some  what  similar  results  have  been 
obtained  for  PBS  9501  (Wiegand  and  Reddingius,  2005a). 

In  the  higher  pressure  range  which  is  about  7.0  MPa  to 
138  MPa  for  EDC37  the  slope  at  larger  strains  is  positive, 
i.e.  work  hardening  is  observed  and  a  maximum  in  the 
stress-strain  curve  is  not  observed  (not  shown).  In  this 
pressure  range  the  initial  slope,  which  gives  the  the 
modulus,  the  yield  strength  and  the  work  hardening  slope 
all  increased  with  increasing  pressure  but  are  much  less 
sensitive  to  pressure  than  the  modulus,  the  compressive 
strength  and  the  work  softening  slope  in  the  low  pressure 
range.  Somewhat  similar  results  have  been  obtained  for 
PBS9501  (Wiegand  and  Reddingius,  2005))  and  LX- 14 
(Wiegand  and  Reddingius,  unpublished  results).  As  noted 
above  the  higher  pressure  range  will  be  discussed  in  detail 
elsewhere  (Wiegand  et  al,  to  be  published). 

In  Figure  2  the  compressive  strength,  the  maximum 
compressive  stress,  is  given  versus  confining  pressure  in 
the  low  pressure  range  for  EDC37,  and  a  straight  line  has 
been  fitted  to  the  data  points.  The  results  indicate  a  linear 
increase  with  increasing  pressure  with  a  slope  close  to  two. 
Similar  results  were  obtained  at  other 


EDC37 
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Figure  2.  Compressive  strength  versus  pressure  for 
pressures  in  the  low  pressure  range. 

strain  rates.  More  limited  results  for  PBS  9501  also 
indicate  a  similar  increase  of  the  compressive  strength  with 
increasing  pressure  (Wiegand  and  Reddingius,  2005).  The 
results  of  Figures  1  and  2  indicate  that  crack  processes  are 
inhibited  by  confining  pressure,  i.e.,  the  stress  required  at 
any  strain  to  cause  crack  damage  increases  with  increasing 
confining  pressure. 

As  noted  above  the  yield  strength  for  EDC37  also 
increases  with  increasing  pressure  in  the  higher  pressure 
range  (not  shown)  .  However,  the  slope  in  the  higher 
pressure  range  is  about  1/40  of  the  slope  of  the 


compressive  strength  versus  pressure  as  given  in  Figure2 
for  the  low  pressure  range. 

After  compression  the  samples  are  barreled,  i.e.,  the 
diameter  at  the  sample  mid  plane  is  greater  that  the 
diameters  at  the  sample  ends.  However,  in  all  cases 
there  is  an  increase  in  the  sample  average  diameter  and  a 
decrease  in  the  average  sample  length  due  to 
deformation. 

The  dimensions  and  densities  of  all  samples  were 
measured  at  atmospheric  pressure  before  and  after 
compression.  The  dimensions  were  determined  by 
measurements  at  several  locations  on  the  sample  and 
averaging.  As  noted  above  the  lengths  of  all  samples 
after  compression  are  decreased  relative  to  values  before 
compression  so  that  all  have  permanent  negative 
(compressive)  axial  strains.  In  addition,  the  average 
diameters  of  all  samples  after  compression  are  increased 
relative  to  values  before  compression  so  that  all  have 
positive  permanent  radial  strains.  The  volumes  and 
densities  of  most  samples  were  determined  by  weighing 
in  air  and  in  purified  water  and  by  use  of  the  water 
density  at  the  temperature  of  measurement.  The  volumes 
before  and  after  compression  indicate  small  volume 
expansions  due  to  deformation.  In  addition  the  density 
measurements  indicate  small  fractional  density  decreases 
with  deformation  ranging  from  hundredths  of  a  percent 
up  to  about  3.9%. 


4.  DISCUSSION 

The  results  in  the  low  pressure  range  can  be 
understood  on  the  basis  of  a  friction  model  developed  by 
Dienes  (1983)  and  Kuo  and  Dienes  (2002).  Other 
authors  have  also  considered  the  effect  of  friction  on  the 
mechanical  properties  of  materials  (McClintock  and 
Walsh,  1962)  (Mavko,  1979)  (Comninou  and  Dundurs, 
1980)  (Tian-Hu  and  Weishen,  1993)  (Egami  and 
Kitaoka,  1996)  (Kovtunenko,  1998).  Dienes  (1983)  and 
Zuo  and  Dienes  (2002)  have  applied  the  concept  of 
coulomb  friction  forces  between  closed  shear  crack 
surfaces  as  resisting  crack  motion  (Ward  and  Hadley, 
1993).  This  is  illustrated  by  the  sketch  of  Figure  3.  aa  is 
the  applied  compressive  stress,  aap  is  the  applied  stress 
on  a  shear  plane,  p,  whose  normal  makes  an  angle  4)  with 
the  applied  stress  direction,  and  ops  and  apn  are  the 
components  of  aap  parallel  and  perpendicular  to  the  shear 
plane  respectively.  These  latter  stresses  are  given  in 
terms  of  cra  and  <])  as 

(1) 


The  net  shearing  stress  on  the  plane  p  is  then 

CJsnet  —  ^Tps  — P  Gpn  (3) 

where  p  is  the  friction  coefficient.  And  if  a  hydrostatic 
pressure  P  is  applied  to  the  sample  in  addition  to  the 
applied  uniaxial  stress  oa,  equation  (3)  becomes 

asnet  =  C>pS  -p(aan  +  P)  (4a) 

=  aa  Cos  (j)  Sin  4>  -  jo,  (aa  Cos  2  (j)  +  P)  (4b) 

by  using  equations  (1)  and  (2). 

Then  equation  (4b)  can  be  rearranged  to  give 

aa  =  (asnet  +  |J.P)/(Sin  (j)  Cos  (j)  -p  Cos  (j)) 

(5) 


Figure  3.  Schematic  of  stresses  for  the  friction  model. 


If  equation  (5)  is  applied  to  the  maximum  of  the  stress- 
strain  curve,  then  asnct  is  the  shear  stress  at  the  maximum. 
Thus  CTa  increases  linearly  with  P  as  observed  (Figure  2)  if 
CTsnet  is  independent  of  pressure  or  is  linearly  dependent  on 
pressure.  For  the  case  of  asnet  independent  of  pressure  Ga 
increases  linearly  with  pressure  with  a  slope  of 


CTpS  =  CTa  Cos  (j)  Sin  (j) 

2 

(7pn  Ga  CoS  (j) 


(2). 


2 

slope  =  |u/(Sin  tj)  Cos  (j)  -pCos  “  (j)) 


(6) 


(9) 


The  angle  <|>  is  determined  by  the  plane  on  which  the 
shear  stress,  asnet,  is  a  maximum.  By  differentiating 
equation  (4b)  with  respect  to  <])  with  constant  aa,  p  and  P 
and  setting  the  differential  equal  to  zero  the  following  is 
obtained: 

Tan  2  cj)  -  2  jixTancj)  -1=0  (7). 

And  by  taking  the  positive  route 

Tan  tj)  =  jl.1  +  ()ii2  +1)1/2  (8). 

Thus,  by  measuring  the  angle  (|>,  p  can  be  determined 
from  equation  (8).  p  can  also  be  determined  from 
equation  (6)  using  the  slope  of  Figure  2  and  the 
measured  value  of  <|>.  Note  that  the  angle  <|>  is 
independent  of  P.  It  is  assumed  here  that  the  samples  are 
isotropic  and  thus  that  all  crack  orientations  exist. 

The  angle  (|>  can  be  obtained  as  the  angle  that  the 
normal  to  the  fracture  surface  makes  with  the  applied 
stress,  aa,  direction  (see  Figure  3).  However,  in  many 
cases  fracture  was  not  observed  but  white  lines  on  the 
sample  surfaces  were  observed  in  most  cases  when 
fracture  did  not  occur.  (Fracture  was  only  observed  at 
atmospheric  pressure  but  the  white  lines  were  observed 
at  atmospheric  pressure  and  at  some  elevated  pressures 
in  the  low  pressure  range.)  These  white  lines  are  taken 
as  the  precursors  of  the  fracture  surfaces  and  they  make 
approximately  the  same  angle  with  the  applied  stress 
direction  as  the  actual  fracture  surfaces.  Therefore,  the 
angle  (|>  has  been  obtained  from  the  angle  of  the  fracture 
surface  in  those  cases  where  fracture  was  observed  and 
from  the  angle  of  the  white  lines  in  cases  where  they 
were  observed.  The  average  value  is  <j)  =  61.8  degrees. 
By  comparison  this  angle  is  45  degrees  when  p  =  0 
(equation  (8)).  For  plastic  bonded  explosives  this  angle 
has  been  observed  to  be  greater  than  45  degree  in  most  if 
not  all  cases,  thus  indicating  the  general  importance  of 
friction,  p  as  obtained  from  equation  (6)  is  0.58,  and  as 
obtained  from  equation  (8)  is  0.66  and  these  are  listed  in 
Table  2.  The  agreement  of  these  two  values  of  p 
supports  the  assumptions  made  in  obtaining  (6).  p 
obtained  from  equation  (6)  is  expected  to  be  the  dynamic 
coefficient  of  friction  and  the  value  obtained  from 
equation  (8)  is  most  probably  also  the  dynamic 
coefficient. 

Dienes  (1983)  and  Kuo  and  Dienies  (2002)  also  give 
a  relationship  between  the  friction  coefficient,  p, 
Poissons’  ratio,  y  and  the  ratio  of  compressive  to  tensile 
strength  as 

Compressive  Strength)/Tensile  Strength),  = 


[2(2-v)]1/2  [(p2+l)1/2  +  p] 


TABLE  2 

Values  of  the  Friction  Coefficient,  p,  Obtained  by  Using 
Various  Methods 

_ P  Method  of  Determining  u 

0.58  Slope  of  Compressive  Strength  vs.  Pressure 
Curve  and  the  Angle  which  the  Failure  Plane  Normal 
makes  with  the  Direction  of  the  Applied  Stress  - 
Equation  (6) 

0.66  The  Angle  which  the  Failure  Plane  Normal 
makes  with  the  Direction  of  Applied  Stress  -  Equation  (8) 

>0.62  The  Ratio  of  Compressive  Strength  to  Tensile 
Strength  -  Equation  (10) 

In  equation  (9)  the  compressive  strength  and  the  tensile 
strength  are  the  threshold  stresses  required  to  initiate  rapid 
unstable  crack  growth  in  compression  and  tension 
respectively  (Dienes,  1983)  (Rice,  1984)  (Dienes,  1984). 
In  tension  the  stress  increases  in  an  approximately  linear 
manner  with  increasing  strain  to  a  maximum  and  then  the 
sample  fails  in  a  brittle  fashion  by  fracture  (not  shown) 
(Ellis,  et  al,  2005)  (Wiegand,  unpublished  results).  In 
addition  there  is  very  little  acoustic  emission  with 
increasing  tensile  stress  until  the  maximum  stress  is 
reached  (Ellis  et  al,  2005).  This  acoustic  emission  is 
interpreted  as  due  to  elastic  waves  generated  primarily  or 
at  least  in  part  by  crack  processes  (Ellis  et  al,  2005).  Thus, 
the  results  suggest  that  there  is  minimal  crack  activity  until 
the  maximum  tensile  stress  is  achieved  and  that  this 
maximum  stress  is  the  stress  necessary  to  initiate  rapid 
unstable  crack  growth.  In  contrast,  in  compression 
acoustic  emission  has  been  observed  throughout  the 
increasing  stress  part  of  the  stress  strain  curve,  through  the 
maximum  and  into  the  work  softening  part  of  the  curve 
(Ellis  et  al  2005),  but  fracture  indicating  unstable  rapid 
crack  growth  has  not  been  observed  (See  the  stress-strain 
curve  at  atmospheric  pressure  in  Figure  1).  It  is  suggested 
that  only  slow  crack  growth  (Charles,  1958)  occurs  during 
the  compressive  stress-strain  curve  and  that  the  observed 
acoustic  emission  is  primarily  or  at  least  in  part  due  to  this 
slow  crack  growth.  It  is  further  suggested  that  the  damage 
introduced  by  this  slow  crack  growth  so  weakens  the 
sample  that  the  stress  required  for  unstable  rapid  crack 
growth  in  compression  is  not  attained.  In  some  plastic 
bonded  explosives,  e.g.,  PAX  2A  rapid  unstable  crack 
growth  is  observed  in  compression  at  low  temperature 
(Wiegand,  1999.)  However,  at  25  C  the  tensile  and 
compressive  stress-strain  curves  for  PAX  2A  are  very 


similar  to  those  of  EDC37.  Slow  crack  growth  is  not 
expected  at  low  temperatures. 

Slow  crack  growth  is  thermally  activated  and  stress 
assisted  (Charles,  1958)  and  does  not  have  an  explicit 
threshold  stress.  In  contrast  rapid  unstable  crack  growth 
does  have  a  threshold  stress.  Dienes  (1998)  and  Dienes 
and  Reilly  (1998)  used  crack  growth  models  without 
thresholds  to  obtain  fits  to  Wiegand’s  (unpublished 
results)  uniaxial  compressive  stress-strain  data  for  PBX 
9501.  In  addition  Dienes  (1998)  was  able  to  obtain  a  fit 
to  the  strain  rate  dependence  of  Wiegand’s  results  only 
by  introducing  slow  crack  growth  into  his  model.  And 
modelers  at  AWE,  Aldermaston  have  found  a  very  low 
stress  threshold  for  crack  growth  in  EDC37.  All  of  these 
results  support  the  hypothesis  that  primarily  slow  crack 
growth  and  not  rapid  unstable  crack  growth  occurs  in 
compression  in  these  composites  for  the  conditions  of 
interest  here.  The  maximum  stress  observed  in 
compression,  the  compressive  strength,  is  then  less  that 
the  stress  required  for  rapid  unstable  crack  growth  and 
this  latter  stress  is  not  observed  for  the  reason  given 
above.  Returning  then  to  equation  (9)  a  calculation  of 
the  friction  coefficient  using  the  compressive  strength 
and  the  tensile  strength  will  yield  a  coefficient  which  is 
less  than  the  true  value  because  the  compressive  strength 
is  less  than  the  stress  required  for  rapid  unstable  crack 
growth.  Poisson’s  ratio  has  been  found  to  be  0.42  for 
EDC37  in  the  low  pressure  range  using  neoprene  shrouds 
and  the  system  described  above  in  the  experimental 
section.  The  ratio  of  compressive  to  tensile  strengths  is 
3.2  for  the  conditions  of  this  work  (Ellis,  unpublished 
results).  Then  |i  obtained  from  equation  (9)  is  0.62  and 
this  is  listed  in  Table  3  as  a  lower  limit.  The  agreement 
between  two  of  the  values  of  p,  Table  2,  as  obtained 
from  equations  (6),  and  (8)  and  the  minimum  value  as 
obtained  from  equation  (9)  provide  significant  support 
for  the  role  of  friction  in  determining  the  compressive 
stresses.  Values  of  p  for  three  other  composite  plastic 
bonded  explosives  obtained  using  equation  (9)  are 
somewhat  close  to  the  values  for  EDC37. 

It  is  to  be  noted  that  at  high  strain  rates  the  time 
available  for  slow  crack  growth  for  a  given  range  of 
strain  will  be  much  less  than  the  time  needed  for  the 
same  range  of  strains  at  low  strain  rates.  Therefore,  the 
damage  introduced  by  slow  crack  growth  will  be  less  at 
the  higher  strain  rate.  Thus  the  sample  will  support 
larger  stresses  at  the  higher  rate  and  so  the  stress  may 
increase  to  that  required  for  rapid  unstable  crack  growth 
and  fracture. 

Because  there  is  friction  between  crack  surfaces 
associated  with  crack  motion,  there  will  be  heating  of  the 
crack  surfaces.  This  raises  the  possibility  of  ignition  of 
reaction  at  crack  surfaces  because  of  this  heating.  Kuo 
and  Dienes  (2002)  have  presented  evidence  for  reaction 


at  cracks  surfaces  taken  from  the  work  of  Howe  et  al 
(Howe  et  al,  1985).  These  latter  investigators  impacted 
artillery  shells  containing  TNT  and  afterwards  polished 
TNT  surfaces  to  reveal  cracks.  In  particular,  several 
cracks  exhibited  evidence  of  reaction  along  the  crack 
length,  thus  giving  support  to  the  hypothesis  of  ignition 
due  to  frictional  heating  at  the  crack  surfaces.  In  addition. 
Dienes  et  al  (2006)  have  calculated  particle  velocity  versus 
time  curves  which  are  very  similar  to  the  observed  curves 
of  Mulford  et  al,  (1993)  for  PBX  9501  subjected  to 
multiple  shocks.  Dienes  et  al  used  a  model  in  which 
frictional  heating  in  shear  cracks  raises  the  local 
temperature  and  pressure  sufficiently  to  initiate  reaction. 
Both  of  these  works  indicate  the  very  practical  importance 
of  friction  between  the  surfaces  of  cracks  in  explosives. 

5.  SUMMARY  AND  CONCLUSIONS 

A  linear  increase  of  the  compressive  strength  with 
increasing  pressure  is  attributed  to  the  effect  of  coulomb 
friction  between  the  surfaces  of  closed  shear  cracks.  The 
slope  obtained  from  this  linear  relationship  is  proportional 
to  the  friction  coefficient  and  is  also  dependent  on  the 
failure  angle,  the  angle  which  the  failure  plane,  the  plane 
of  maximum  shear  stress,  makes  with  the  loading 
direction.  In  addition,  Dienes  (1983)  and  Kuo  and  Dienes 
(2002)  predict  a  relationship  between  the  friction 
coefficient  and  only  the  failure  angle.  The  friction 
coefficient  obtained  from  the  slope  of  the  strength  versus 
pressure  curve  and  the  failure  angle  is  in  substantial 
agreement  with  the  coefficient  obtained  only  from  the 
failure  angle.  These  authors  also  predict  a  relationship 
between  the  friction  coefficient  and  the  ratio  of 
compressive  strength  to  tensile  strength,  at  atmospheric 
pressure.  A  lower  limit  on  the  value  of  the  friction 
coefficient  has  been  obtained  from  this  relationship  which 
is  in  agreement  with  the  other  values  noted  immediately 
above.  This  agreement  of  friction  coefficients  obtained 
from  the  three  different  relationships  is  taken  as  strong 
support  for  the  role  of  friction  in  the  mechanical  properties 
of  this  composite. 
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